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FRACTAL IMAGE DATA AND IMAGE GENERATOR 
FIELD OF THE INVENTION 

The present invention relates to fractals, and in particular both to a new class of fractals 
5 and also to a fractal generator able to iteratively generate image data and a number of such 
fractals. 

BACKGROUND 

Since the 1980's, fractals have been used to describe a wide variety of physical systems. 
Fractals can be generally classed as deterministic, if they can be generated by a 

10 deterministic process, or random, if the generation process uses one or more probability 
distributions to select iterated function systems. One of the best known applications of 
fractals is the simulation of natural formations occurring in the physical world, such as 
mountains, clouds, and waves. One of the difficulties of existing fractals and procedures 
for their generation is that deterministic fractals tend to be too rigidly deterministic, and 

15 random fractals too random, when these are compared with natural processes and forms 
occurring in the physical world. Deterministic fractals, while initially eye-catching, tend to 
lose the viewer's interest because they appear too repetitive. On the other hand, random 
fractals can also appear somewhat boring because the eye finds it hard to discern patterns, 
similar to the way that wispy white clouds can appear boring. These weaknesses are one of 

20 the factors that have limited the use of fractals in the content industries, in the areas of 
special effects in digital animation, computer games, and advertising, for example. 

It is desired to provide at least a useful alternative, and in particular to provide image data 
representing at least a nascent fractal and a generator and process for generating the image 
25 data that alleviate one or more of the above difficulties. 
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SUMMARY OF THE INVENTION 

In accordance with the present invention, there is provided a V-variable fractal. 

5 In particular, the present invention provides a V-variable fractal represented by fractal 
image data, where V is an integer greater than one and represents the number of 
constituent images available for iterative combination to generate the fractal. The images 
are combined in a random manner. 

10 The present invention also provides image data representing a variable number n of 
constituent images iteratively transformed and combined in a random manner to generate 
said image data, with 1 < n < V. 

The present invention also provides a fractal generation process, including: 
1 5 (i) randomly selecting images from a set of input images; 

(ii) selecting transformation functions from a set of transformation functions; 

(iii) generating transformed images by applying the selected transformation 
functions to the selected images; and 

(iv) generating an output image by combining the transformed images. 

20 

The present invention also provides a fractal generation process, including randomly 
selecting from a set of input images, transforming the selected images, and combining the 
transformed images to generate a set of output images, and iterative repetition of these 
steps using the set of output images of each iteration as the set of input images for the next 
25 iteration. 



WO 2004/081876 



PCT/AU2004/000312 



-3- 

The present invention also provides a fractal generator, including: 

an image selector for selecting M images from V input images; 
a function selector for selecting a set of M transformation functions; 
at least one image transformer for respectively applying the selected transformation 
5 functions to the selected input images; and 

a compositor for composing an output image from the images output by said at 
least one image transformer. 

The present invention also provides a fractal generation system, including an image 
10 selector for selecting images from a set of input images, and an image transformer for 
transforming the selected images to generate a set of output images, said system being 
adapted to provide said set of output images as the set of input images to iteratively 
generate fractal image data. 

15 The present invention also provides fractal image data representing a combination of two 
or more constituent first images, each of said first images representing a random 
transformed combination of two or more constituent second images, each of said second 
images representing a random transformed combination of two or more constituent third 
images, each of said third images representing a random transformed combination of two 

20 or more constituent fourth images, wherein each transformation includes at least one of 
translation and rotation. 

The present invention also provides image data decomposable into at least four successive 
levels, wherein each level is composed of smaller data sets which are affine 
25 transformations of V basic sets. 

The present invention also provides image data representing iterative transformation and 
combination of at least two images selected from a set of V > 1 input images, wherein 
image data generated at each iteration represents a combination of at least two smaller 
30 images, wherein each of said at least two smaller images represents an affine or projective 
transformation of image data generated at the previous iteration. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred embodiments of the present invention are hereinafter described, by way- of 
example only, with reference to the accompanying drawings, wherein: 

Figure 1 is a block diagram of a preferred embodiment of a fractal generation 
system; 

Figure 2 is a block diagram of a fractal generator of the fractal generation system; 
Figure 3 is a block diagram of an image transformer of the fractal generator; 
Figure 4 is a flow diagram of a fractal generation process executed by the fractal 
generation system; 

Figure 5 is a schematic diagram illustrating transformation of an input image by a 
transformation module of the image transformer; 

Figures 6 to 10 are screenshots of images generated by the system using a first set 
of input parameters; 

Figures 11 to 15 are screenshots of images generated by the system using a second 
set of input parameters; and 

Figures 16 to 23 are screenshots of images generated by the system using a third set 

of input parameters. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

As shown in Figure 1, a fractal generation system includes a fractal generator 100. As 
shown in Figure 2, the fractal generator 100 includes an image selector 202, an iterated 
function system (IFS) selector 204, and an image transformer 206. The fractal generation 
system executes a fractal generation process that generates output image data 102 from 
input image data 104, iterated function systems 106, and selection probabilities 108. The 
input image data 104 constitutes a set of input images 104, and the output image data 102 
constitutes two or more output images 102 that can be considered as representing two or 
more fractals or nascent fractals, as described below. 



WO 2004/081876 



PCT/AU2004/000312 



-5- 

In the described embodiment, the fractal generation system is a standard computer system 
such as an Intel™ x86-based personal computer including a Pentium™ processor 110, 
random access memory (RAM) 112, and non-volatile (e.g., magnetic disk) storage 114, 
and the fractal generator is implemented by software modules stored on the non-volatile 
5 storage 114 of the computer and executed by the processor 110, However, it will be 
apparent that at least parts of the fractal generator can be alternatively implemented by 
dedicated hardware components, such as application-specific integrated circuits (ASICs). 

The fractal generation system is described as generating output image data 102 from input 
10 image data 104. Although it is not necessary that the input and output data relate to images 
as such, i.e., the system processes these as arbitrary numeric data, it is expected that the 
system will be used to process input data and generate output data that in each case will be 
visualised as one or more images on a display or hardcopy. Accordingly, the input image 
data 104 is hereinafter referred to as input images 104, and the output image data 102 is 
1 5 hereinafter referred to as output images 1 02. 

A user of the fractal generation system provides to the system a number V > 1 (i.e., at least 
two) input images 104, together with N > 1 iterated function systems F = {F 1 , F N } 106 
and N selection probabilities P = {P 1 , . . P N } 108. Each of the N iterated function systems 

20 106 comprises a set of M transformation functions. Each transformation function defines a 
contractive affine transformation that can be applied to an image, and is represented by the 
six parameters r, s, x, y,0 and % as described below. Alternatively, other types of 
transformations can be used, each represented by a corresponding set of parameters. As 
shown in Figure 5, an input image 502 (in this case having horizontal and vertical 

25 dimensions of unity) is transformed by first scaling the horizontal and vertical dimensions 
of the input image 502 by the factors r and s, respectively, with -1 < r < 1 and -1 < s <1, 
followed by translation along the x and y axes by the amounts specified by the x and y 
parameters, respectively, followed by a geometrical distortion whereby horizontal lines are 
rotated by the angle 0, and vertical lines are rotated by the angle <p, with -180 < 6 < 180 

30 and - 180 < <p < 180 . In Figure 5, the angles 6 and q> are nearly equal, and the output 
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image 504 appears to retain its rectangular shape. However, when these two parameters 
are substantially unequal, a rectangular input image 502 is distorted to provide a generally 
parallelogram or diamond-like shape. 

5 Because each transformation function is represented by the six transformation parameters 
described above, each iterated function system (IFS) 106, being a set of M transformation 
functions, can be represented as a parameter table 506 with M rows, each row representing 
a distinct transformation function. For example, an IFS 106 for M=3 can be represented as 
a 3x6 matrix or table, as follows: 

10 



ri 


Si 


Xl 


yi 


9i 


<pi 




S2 


x 2 


y2 


e 2 


q>2 


T3 


S3 


x 3 


ys 


e 3 


<P3 



Thus the iterated function systems 106 are provided as a set of N> 1 parameter tables, 
each having M > 1 rows. The N selection probabilities 108 are used to select a particular 
IFS (i.e., one of the N parameter tables) to apply to selected images, as described below. 

15 

As shown in Figure 2, the fractal generator 100 operates on a set of V input buffers 208 to 
212, and a set of V output buffers 214 to 218 allocated from the computer RAM 112. 
Although Figure 2 shows V=3 of each kind of buffer, there can be in general an arbitrary 
number V > 1 of each buffer type, subject only to system constraints. The number of input 

20 images 104 provided to the system is less than or equal to the number V of input and 
output buffers, as the same image can be allocated to more than one buffer. The input 
images 104 are copied into the input buffers 208 to 212, and the fractal generator 100 then 
executes a fractal generation process, as shown in Figure 4. The fractal generator 100 
generates output images into each of the output buffers 214to218ina sequential manner, 

25 starting with the first output buffer 214, and ending with the last output buffer 218. When 
the outermost loop of the fractal generation process is executed, all the output buffers 214 
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to 218 have been updated with new images, and an iteration of the process is said to have 
completed. 

The fractal generation process begins at step 402, when the first output buffer 214 is 
5 selected to receive the first output image. At step 404, the image selector 202 randomly 
selects M of the V input images stored in the input buffers 208 to 212, with the same image 
allowed to be selected more than once. Accordingly, each time a selection is made, the 
probability of selecting the input image from a particular buffer is equal to 1 /V, although it 
will be apparent that alternative selection methods can be used. The M selected images are 

10 provided to the image transformer 206. At step 406, the IFS selector 204 randomly selects 
one of the N IFSs 106 using the respective N selection probabilities 108. However, the 
random selection can be a biased quasi-random selection if desired. For example, the N 
selection probabilities 108 can be provided as respective weights used to divide the 
numeric interval between 0 and 1 into N contiguous subdivisions, where .the size of each 

15 subdivision is proportional to the weight for that subdivision. A random number between 
0 and 1 can then be generated to select one of the N IFSs 106, according to which 
subdivision the random number falls within. 

As shown in Figure 3, the image transformer 206 includes M transformation modules 302, 
20 304, and a superimpose module 306, shown for the example of M=2. At step 408, the 
transformation modules 302, 304 transform the selected input images using the selected 
IFS. The selected input images are transformed using respective transformation functions 
from the selected IFS. That is, the first transformation module 302 applies a first 
transformation function defined by the first row 312 of the selected parameter table 316 to 
25 the first selected input image 308. Similarly, the second transformation module 304 applies 
a second transformation function defined by the second row 314 of the selected parameter 
table 316 to the second selected input image 310. 

Each of the transformation modules 302, 304 generates a transformed version of its input 
30 image, using the corresponding transformation function, as described above. It will be 
appreciated that, although distinct transformation modules 302, 304 are shown in Figure 3, 
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the M selected input images could alternatively be processed sequentially by a single 
transformation module. However, it may be particularly advantageous to provide distinct 
execution instances of each of the transformation modules 302, 304 if the system includes 
more than one processor or if the transformation modules 302, 304 are implemented as 
5 dedicated hardware components. 

In the preferred embodiment, each of the input images 104 provided to the system by a 
user is a rectangular image with dimensions m x n utilising a rectangular co-ordinate 
system with origin 0,0 corresponding to one of the corners of the rectangular input image. 

10 Each of the transformed images generated by the transformation modules 302, 304 is 
maintained as an m x n image corresponding to the same spatial co-ordinates as the 
original input images 104. Accordingly, a given transformed image will generally include 
a number of transformed versions of previous images positioned within an otherwise 
empty mxn space within the particular transformed image. Any transformed images or 

15 parts of transformed images that fall outside the mxn co-ordinate space of a transformed 
image are cropped and are not included within the transformed image. Because each 
transformed image corresponds to the same m x n region of coordinate space, transformed 
images can be combined by superposition. However, it will be apparent that alternative 
embodiments can be devised such that a transformed image is not generated as an m x n 

20 rectangular image within the m x n coordinate space defined by input images 104, and the 
combining of images may involve only partial overlap of the constituent images, or may 
not involve any overlap at all. 

At step 410, the resulting transformed images are processed by the superimpose module 
25 306. This involves superimposing or overlaying the transformed images to generate a 
single output image 318. The superimposition is achieved by assigning to each pixel of the 
output image 318 a colour value (e.g., comprising a triplet of red, green, and blue (RGB) 
colour components) equal to the average of the colour values of the corresponding pixels 
of the transformed images. The resulting superimposed output image 318 is stored in the 
30 current output buffer, in this case the first output buffer 214. 
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It will be apparent that the transformed images can be combined in a variety of alternative 
ways. For example, the images can alternatively be combined by putting one down on top 
of the preceding one, one after another, so that only the last colour value written to each 
pixel is kept. Alternatively, a depth can be associated with each pixel in each image, for 
example, by adding another dimension to each image and increasing the parameter set 
from 2D to 3D, and then, when combining images, the frontmost pixel from the images 
being combined is the one that is used, as is done in 3D computer graphics, where a Z- 
buffer is used to select the frontmost pixels. Or alternatively, the exclusive or (XOR) of the 
first image and the second image can be used. 

At step 412, a check is performed to determine whether the current output buffer is the last 
output buffer 218. If not, then the next output buffer, in this case the second output buffer 
216, is selected at step 414, and the process returns to step 404 to randomly select another 
15 M images from the V input buffers 208 to 212, and subsequently generate another 
superimposed output image 318 for storage in the second output buffer 216 by executing 
steps 406 to 410. These steps are repeated until all the V output buffers 214 to 218 have 
been used to store respective new superimposed output images 318. Once the last output 
buffer 218 has stored an image, the test at step 412 succeeds, and an iteration of the 
20 process is complete. Rather than selecting another output buffer, the contents of the output 
buffers 214 to 218 are then used as input images to the next iteration of the process. This 
can be achieved by physically copying the contents of the output buffers 214 to 218 into 
the respective input buffers 208 to 212. However, it is more efficient simply to exchange, 
at step 416, the roles of the two sets of buffers by exchanging or swapping the 
25 corresponding buffer pointers into the random access memory 112. In any case, after the 
superimposed output images have been provided in respective input buffers 208 to 212 and 
the output buffers 214 to 218 have been cleared, the fractal generation process begins the 
next iteration by selecting the first output buffer 214 at step 402, and repeating the steps 
described above. As many iterations as desired can be performed. 



30 
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The contents of one or more of the output buffers 214 to 218 can be displayed on a display 
device associated with the fractal generation system for viewing by the user. Individual 
images selected by the user or images generated by a user-specified number of iterations of 
the fractal generation process can be stored as final output images 102 on the non-volatile 
5 storage medium 114 for subsequent use. After a number of iterations, the output images 
102 are a set of V-variable fractal images derived from the preceding input images. It will 
be appreciated by those skilled in the art that the although the image data generated by the 
fractal generation system is referred to as fractal image data, it may only represent a 
nascent fractal if the finite number of iterations of the fractal generation process are not 
1 0 sufficient to produce a fractal. 

With successive iterations, the output images 318 become less dependent on the initial 
input images 104, and eventually become characteristic of the input parameter tables 106 
and their selection probabilities 108. The images 318 generated by the system then form a 

15 diverse sample of an infinite collection of all possible images associated with the initial 
parameter tables 106 and probabilities 108. In particular, because the transformations are 
contractive, the original input images 104 for the first iteration will eventually be 
contracted to the size of a single pixel when the entire output image is viewed, given a 
sufficient number of iterations. Thus the output images become independent of the 

20 original input images 104. Accordingly, the input images 104 can be alternatively 
provided or generated by the system itself, rather than being provided by the user as input, 
and the value of V can be provided as input to the system. By adjusting the parameter 
tables 106, the appearance and form of the classes of images generated by the system can 
be modified. For example, one class of images might resemble a group of tea trees, 

25 another a collection of faces, and another might resemble (and be vised to represent) 
possible stock market simulations. As a result, the possible applications of the fractal 
generation system and process are very broad and diverse. 

For example, Figures 6 to 10 are screenshots of images generated by the fractal generation 
30 system after 4, 5,6, 11 and 44 iterations, respectively, using V=2 and the following input 
parameter tables: 
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IFS 1: 



transformation 1: 
transformation 2: 



a 


b 


c 


d 


e 


f 


0.5 


-0.375 


0.3125 


0.5 


0.375 


.1875 


0.5 


0.375 


0.1875 


-0.5 


0.375 


.6875 



5 IFS 2: 





a 


b 


c 


d 


e 


f 


transformation 1: 


0.5 


-0.375 


0.3125 


-0.5 


-0.375 


.8125 


transformation 2: 


0.5 


0.375 


0.1875 


0.5 


-0.375 


.3125 



In this example, each transformation is a contractive affine transformation defined by 
( x> y) _> (ax+by+Q^ dx+ey+f), where the parameters a to f are as indicated in the parameter 
tables. 

As a second example, Figures 11 to 15 are screenshots of fractal images generated after 1, 
3, 6, 8, and 22 iterations, respectively, using V=2 and the following parameter tables: 



IFS 1: 



a 


b 


c 


d 


e 


f 


g 


h 


I 


.151827 


.034074 


-.215212 


-.123706 


-.246197 


.054262 


-.212523 


.059302 


-.234997 


.011248 


.011321 


-.000425 


-.001814 


-.026793 


.062591 


.014512 


-.057923 


.086631 


3 2: 

a 


b 


c 


d 


e 


f 


g 


h 


I 


.017380 


-.035039 


.046965 


.029178 


.006695 


.025663 


.021377 


-.012129 


.048469 


.144191 


.004070 


-.183496 


-.099432 


-.243409 


.043246 


-.151388 


.001065 


-.199968 
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In this example, each transformation is a contractive projective transformation defined by 
(x,y) -» ( (ax+by+c) / (gx+hy+I), (dx+ey+f) ) / (gx+hy+I) ) , where the parameters a to I 
are as indicated in the parameter tables. 

5 As a third example, Figures 16 to 23 are screenshots of images generated after 1, 2, 3, 4, 6, 
8, 22, and 40 iterations, respectively, using the same projective transformation with V=2 
and the following parameter tables: 



IFS 1: 



a 


b 


c 


d 


e 


f 


S 


h 


I 


.001231 


-.335970 


.002318 


-.164124 


-.028044 


-.125233 


.003548 


-.033609 


-.309261 


.002575 


.039210 


.007817 


-.045553 


-.005329 


.056137 


-.001851 


-.027965 


.126001 


S2: 


a 


b 


c 


d 


e 


f 


g 


h 


I 


-.006211 


-.197079 


.208860 


.087110 


.013677 


.011083 


-.008298 


.075200 


.205408 


.003412 


.306195 


-.303822 


.141253 


.023639 


-.286545 


.001890 


.028936 


-.312419 



The image data or images generated by the system represent a new type or class of fractals, 
referred to as V-variable fractals. As can be seen from Figures 6 to 23, they can be 
15 considered to be random fractals because they exhibit a degree of randomness, but this is 
intermediate between deterministic fractals and prior art random fractals generated by 
other methods. As such, they are able to represent natural forms more accurately than prior 
art fractals, and yet can be generated more readily than prior art random fractals. 

20 An image representing a V-variable fractal can be characterised as follows. Each fractal 
image is a combination of at least two smaller images. There is a set of V basic images 
such that each of the smaller images represents a contractive affine transformation of a 
basic image from the set of V basic images. The V basic images and the contractive affine 
transformations provide the first level of decomposition of the image representing the V- 

25 variable fractal. Each of the two or more smaller images from the first level of 
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decomposition is itself a combination of at .least two smaller again images and there is a 
second set of V basic images such that each of the smaller again images represents a 
contractive affine transformation of a basic image from the second set of basic images. 
The second set of V basic images and the corresponding contractive affine transformations 
5 provide a second level of decomposition of the image representing the V-variable fractal. 
These decompositions can be continued for at least four levels of decomposition. However, 
if every set of basic images contains just one distinct basic image, then the representing 
image is a deterministic or random fractal of a previously known type, and as such is not 
included in the new class of fractals. 

10 

It will be apparent to those skilled in the art that the affine transformations in the above 
description can be replaced by projective transformations or certain other classes of 
transformations. These classes have the property that when any two transformations from 
the class are applied one after another, a new transformation in the same class is obtained. 

15 

Figures 6 to 10, 13 to 15, and 19 to 23 are examples of V-variable fractal images where 
V=2. Each of these images was generated by four or more iterations of the fractal 
generation process. At each scale there are at most V non-equivalent sets. These sets can 
be combined partially or wholly one on another. The sets typically change from 
20 magnification level to magnification level. While V=2 in Figures 6 to 23, larger values of 
V can be used. 

The properties of V-variable fractals are described in more detail in M. Barnsley, 
J. Hutchinson and 0. Stenflo, A Fractal Valued Random Iteration Algorithm and Fractal 
25 Hierarchy, available at http://arxiv.org/abs/math.PR/03 12187 : and in M. Barnsley, 
J.Hutchinson and 6. Stenflo, V-variable fractals and superfractals at 
http://arxiv.org/abs/math.PR/03 123 14 . 

The fractal generation system has been described above in terms of contractive affine 
30 geometrical transformation and superimposition of two-dimensional images. However, it 
will be apparent that the iterated function systems 106 can include other function systems 
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that can be used to modify the input image data 104. For example, the iterated function 
systems 102 can include contractive projective transformations, or affine transformations 
that are contractive on average, where the average is determined by the probabilities 108 
associated with the iterated function systems 106. Because the system generates output 
5 images by combining multiple, transformed copies of input images, contractive 
transformations are preferably used in order to limit the output image size when translated 
and/or rotated images are combined, and this is preferably limited to the size of the input 
images. However, the transformations can alternatively be contractive on average, as 
described above. Similarly, it will be apparent that the transformed images can be 
10 combined by methods other than superposition, as described above. Moreover, the fractal 
generation process can be applied to image data representing three-dimensional images to 
generate a new class of three-dimensional fractals, or indeed to input data of even higher 
dimensionality. The image data need not represent an existing image or be used to 
generate an image, but can be used for any desired purpose. 

15 

In an alternative embodiment, the transformation modules 302, 304 can modify the colour 
and/or brightness of input images, in addition to the geometrical transformations described 
above. These additional modifications are defined by additional parameters included in the 
parameter tables 106. When the modifications of colour and/or brightness are performed in 

20 a contractive manner, the output images generated by the system, after a small number of 
iterations, depend only on the input parameter tables 106 and their selection 
probabilities 108. If the brightness of an image is represented by a number z, then an 
example of a contractive transformation of brightness z is new_z= p*z + q where |p|<l and 
q is a constant. Contractive transformations of colour are achieved by three contractive 

25 transformations of brightness, one for each of the red, blue and green colour components. 

In a further embodiment, the parameter tables are generated graphically. For example, an 
affine transformation can be characterised by how it transforms a given triangle into 
another triangle. For example, a user can use a pointing device (such as a mouse) to drag 
30 the vertices of a given triangle to new positions and thereby define an affine 
transformation. 
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In yet a further alternative embodiment, the parameter tables 106 and selection 
probabilities 108 are themselves generated from predetermined probability distributions, 
such as normal distributions, using a random number generator. 

5 

It may be appreciated that the new class of V-variable fractals, realised as images or in any 
digital form, may also be generated by processes other than those described above. For 
example, those skilled in the art will realise that for the generation processes described 
herein, which correspond to forward methods, there are also corresponding backward 
10 methods, although such methods will be less efficient for generating V-variable fractals. 

Many modifications will be apparent to those skilled in the art without departing from the 
scope of the present invention as herein described with reference to the accompanying 
drawings. 

15 



